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The approach pictures taken by the Viking 1 and Viking 2 spacecrafts two 
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Mars data analysis project in order to search for new satellites within 
the orbit of Phobos. To accomplish this task, a complete search proce- 
dure and analysis strategy were formulated, developed and executed using 
the substantial image processing capabilities of the Image Processing 
Laboiatory at the Jet Propulsion Laboratory. The development of these 
new search capabilities should prove to be valuable to NASA in processing 
of image data obtained from other spacecraft missions. The result of 
applying the search procedures to the Viking approach pictures was as 
follows: no new satellites of comparable size (-20 km) and brightness 

to Phobos or Demios were detected within the orbit of Phobos. 
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I . INTRODUCTION 


The objective of this Mars data analysis project has been to search 
the Viking 1 and Viking 2 approach pictures for new satellites. In 
executing this objective, the project has utllifsed the substantial image 
processing capabilities of the Image Processing Laboratory (IPL) at the 
Jet Propulsion Laboratory, the support for which was supplied separately 
by NASA Headquarters. 

The original Interest In pursuing this project csto from data 
obtained on April 8, 1976 by McCook (1979) who observed the Mars occul- 
tatlon of the star Epsilon Gemlnlorum. Initial examination of this 
Ingress light curve, recorded at the observatory of Vlllanova Univer- 
sity, showed an intensity drop to total occultation of the light of 
the star for about 0.7 seconds with no apparent change in the background 
light level coming from Mars. This change was associated with some 
occurence in the near-Mars environment and was tentatively interpreted 
as an occultation of the star light by a close Martian satellite with 
perhaps an 8 to 15 km diameter. A later and more careful Interpretation 
of other occultation data (French, Goguen, Duthie, 1978; Llller, et al., 
1978; French and Elliot, 1979) revealed that a similar signature was 
also seen by other observers, but that the signature was most likely 
caused by large-scale structures in the Mars atmosphere of large hori- 
zontal extent (French and Elliot, 1979). 

In order to search the Viking 1 and Viking 2 approach pictures 
for new safellites, a complete search procedure and analysis strategy 
had to be formulated, developed and executed. Since no such capability 
existed at IPL, a considerable effort was undertaken in this project to 
create this needed capability. The search procedure so developed can 
be divided into three parts (1) screening of the image data to obtain 
initial satellite candidate objects, (2) navigation processing of the 
candidate object data, and (3) analysis of the candidate object data. 
Each of these parts is described in detail in Section 3. The most dif- 
ficult task was that of developing the image processing techniques for 
the correlation analysis of candidate objects in part 3. This develop- 
ment, however, provided some rather sophisticated software tools and 


capabilities that should prove to be valuable to NASA in processing of 
image data obtained from othe;: spacecraft missions. 

The result of applying the search procedure to the Viking approach 
pictures was as follows: no new satellites of comparable size (~20 km) 

and brightness to Phobos or Deimos were detected within the orbit of 
Phobos. Probing the approach pictures for satellite significantly 
smaller or dimmer than Phobos was not possible because the exposure 
times for the pictures were adjusted to properly Image Mars (an exten- 
ded object) and were not set at their maximum values as would be opti- 
mum for detection of new satellites. 



II. DESCRIPTION OF THE VIKING APPROACH PICTURES 


The Viking 1 and Viking 2 spacecrafts, respectively, were launched 
on 20 August and 9 September in 1975 and initiated their Mars orbital 
Insertion maneuvers on 19 June (l.e., day 171 of the year) and 7 August 
(l.e., day 220 of the year) in 1976. Two days prior to the Viking 1 
orbital insertion maneuvers (i.e., days 169 and 170), nineteen sequences 
of pictures were taken in two or three colors by the two vldicon cam- 
eras. Twc days prior to the Viking 2 orbital insertion maneuvers (i.e., 
days 218 and 219), eighteen sequences of pictures were taken in three 
colors by the two vidicon cameras. These approach plcture£ contained 
all or some fraction of the half -Illuminated disc at Mars. They also 
contained projected views on the image plane of a significant region 
of the surrounding clrcumplanetary space (image-plane displacements 
from the planet center varying from 17 >000 km to 8,000 km for Viking 
1 and from 18,000 km to 7,000 km for Viking 2). One of the Viking 1 
approach pictures, Illustrating the image-plane view of the camera, 
is shown in Figure 1. 

The two known satellites of Mars, Delmos and Phobos, have essentially 
circular orbits that lie nearly in the equatorial plane of the planet 
and have orbital radii of ab^ut 23,000 km and 9,300 km respectively. 

During the time that the approach pictures were taken, both spacecrafts 
were located above south Mars latitudes between 30° and 36° so that their 
cameras imaged the space that was almost always well within the nearly 
equatorial orbit of Delmos and many times well within the nearly equa- 
torial orbit of Phobos (see Figure 1). The approach pictures therefore 
Imaged a region of clrcumplanetary space about Mars suitable to search 
for new satellites near are within the orbit of Phobos. Since the 
orbital period of Phobos is 7.64 hours and the period of an object 
with a circular orbital radius of 4000 km (slightly larger than the 
3393 km equatorial radius of Mars) is 2.13 hours, the approach pictures 
taken mostly at two-hour intervals also have well placed time windows 
in which to detect new satellite within the orbit of Phobos. 

The approach pictures for Viking 1 and Viking 2 are each labeled 
by a unique picture number (PIC. NO.). The picture number is composed 
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of a three-digit number (denoting the day of the year on which the 
picture was taken) followed by a letter (C for Viking 1 and D for 
Viking 2) and finally followed by a sequence number of two digits 
initialized to 01 at the beginning of a new day. The picture numbers 
for each approach picture and some other relevant picture information 
taken from the Viking SEDR (supplementary experimental data records) 
are given in Table 1 and Table 2 for Viking 1 and Viking 2« respectively. 
For Viking 1* two non-siMared approach pictures 168C09 and 168C11 
taken a day earlier than those catalogued in Table 1 aa well aa picture 
169C03 that is listed in Table 1 were mistakenly not Included in our 
study. 

From Table 1 and Table 2, the distance of the spacecrafts from 
Mars during the approach picture exposures can be seen to vary between 
about 700,000 km and 200,000 km. The exposed area on the image plane 
is 1056 pixels vertically (the line value L) and 1182 pixels horizont- 
ally (the sample value S) . The equivalent spatial pixel resolving 
power of the image, the expected size of Phobos and the absolute bright- 
ness of Phobos for these spacecraft distances from Mars are summarized 
in Table 3. A nominal diameter of Phobos of 20 km was adopted in Table 
3 and differs a little from the actual ellipsoid radial dimensions in 
kilometers of (13.5, 10.5, 9.0). The full size of Phobos would vary 
from about one pixel to several pixels, but since it is only about half 
illuminated by the sun, its image signature will appear even smaller. 

The response of the Viking camera to a point light aource is to 
spread the signal from the central pixel to the adjacent pixels so as 
to Increase the image size of the point while reducing its absolute 
brightness per pixel. Phobos, or an object of equivalent brightness, 
should thus produce a distinct several pixel signature on the picture 
image if the exposure time and filter transmission factor (i.e., de- 
termined by the filter color L'«elected, see Table 4) are such to allow 
a sufficient signal to noise level to be captured. From Table 1 and 
Table 2, both the exposure time and filter color can be seen to vary 
from'plcture to picture. The distribution and relative sensitivities 
of these approach pictures are summarized in Table 5, from which it 
may be seen that the Viking 1 pictures are significantly more sensitive 
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than the Viking 2 pictures. For Phoboe, or an equivalently bright 
satellite, the albedo and the photometric function of Its surface 
further reduces the absolute camera response of Table 4 by an amount 
«rhlch Is approxijnately the same for each filter. The ability of the 
cameras to detect Phobos In the approach pictures will be assessed 
and shown in the following section to be possible, unfortunately, 
only for a small subset of the total pictures In Table 1 and Table 2. 
This occurs primarily because the cas^ra exposure values were selected 
to properly Image Mars (an extended object) and were necessarily sig- 
nificantly lower than the maximum exposure times of 2660 msec that 
would be optimum for Identifying new satellites. 


111. SATELLITE SEARCH PROCEDURE AND ANALYSIS 


3.1 In troduction 

In order to utilize the Viking 1 end Viking 2 approach picturee 
listed in Table 1 and Table 2 to search for new satellites* a search 
procedure and analysis had to be fonmilatad and all neceasary iaage 
processing software toola had to be developed. This fomulatlon and 
developnent phase of the project represented a considerable effort. 

It was undertaken at the IPL facility of the Jet Propulsion Laboratory 
(JPL) by Gary Yagi and was supervised by W.H. Smyth of AER. One of 
the direct and positive contributions of this project has been the 
development of a satellite search capability for IPL. This capability 
may prove to be useful in processing of image data from other NASA 
spacecraft missions. 

The satellite search procedure and analysis has been divided into 
three steps: (1) screening of the image data to produce an initial 

candidate object list* (2) first order navigation processing of the 
candidate object data* and (3) analysis of the candidate object data. 

All numerical computations required in this effort were perforiMd on 
the IBM 370/158 computer at IPL and ran under VICAR* the JPL image 
processing language and its associated application programs. In addi- 
tion to appllcati^*n programs under VICAR* many new application programs 
were created* the most diff:l<.ult of which were those needed for the 
prograde and retrograde orbit correlational analysis in step three. 

3.2 Screening of the Image Data 

This first step in the search procedure and analyais is diagrama- 
tlcally illustrated in Figure 2 and will be discussed in detail below. 
The approach pictures are obtained in digital form from the Viking 
Orbiter EDR (experiiMntal data record) data which are stored on compu- 
ter tapes. In order to search for the pixel signature of possible 
satellite candidate objects* a VICAR program was utilized to determine 
the medium background level above which a candidate object could be 
Identified. The background filter chosen for this utilized an 11 x 11 
pixel are., to define the Mdium background at each pixel image location. 



Any object that vet then two or aore adjacent plxela In aite with a 
brightneea of 4 DN or aore above thia background level vaa claaalfled 
as a possible candidate object. A VICAR prograa waa then applied to 
the digital image data to identify and locate the coordinates of the 
centroid of each candidate object in the pixel coordinates (L,S) of 
iaage apace (the geoaetrlcally micorrected iaage plane). The collec- 
tion of these candidate object coordinates for each approach picture 
was then stored on a disc file (see Figure 2). 

Many objecta that are identified by the object locater were caiaera 
blealshea, reseau (uniformly spaced points on the iaage plane) or local 
aaxlauas that occur in the camera system and therefore need to be 
screened from the candidate list. Other objects occurring on the 
disc of Mars as bright spots also need to be so identified and removed 
from the candidate list. Finally transmission errors and other spurious 
signals need to be properly identified and deleted from the candidate 
list. To accomplish this, two methods were applied. A blemish filter 
prograa, which compared known camera blemish patterns with the objects 
in the candidate list, vaa written and applied to the candidate list 
file. A computer printout of the pixel listing of each object that 
succeeded in passing through this blemish filter was then obtained. 
Parallel to this, VICAR programs were used to obtain photographs of all 
approach pictures containing candidate objects with each candidate 
object identified on the photograph. One such photograph for the Vik- 
ing 1 and Viking 2 approach pictures is shown in Figure 3 and Figure 4, 
respectively. Visual examination of the photographs and pixel listings 
was then used to manually screen from the candidate list unsuitable 
objects. Objects rmsalnlng after this screening process define the 
initial candidate object list for the second atage of processing. 

The screening process provided 54 initial candidate objecta for 
the Viking 1 approach pictures and 355 initial candidate objects for 
the Viking 2 approach pictures. The Viking 2 approach pictures contain 
flw>re than six times as many Initial candidate objects as Viking 1. This 
resulted because both the gradient in the absolute background DN-v«lue 
as well as the noise fluctuation about this background value were sig- 
nificantly larger in the Viking 2 approach pictures than those in the 
Viking 1 approach pictures. This difference in the background gradient 
can, for example, be seen by comparing Figure 3 and Figure 4. 


In addition to this ocrooning procooaing, each photograph of tha 
approach pictures that still contained candidate objects was used to 
nanually Mssure the pixel location of the planet center (see Figure 
2). This inforsiation is vital to locate the vector direction of the 
optical axis of each laage which is inaccurately supplied by the cone 
and clock angles in the SEDR. This inaccuracy results from the llait 
cycle of the spacecraft (the uncertainty in the vector orientation of 
the spacecraft that is controlled within a specified tollerance by the 
use of small control jets). 

3.3 First Order Naviaation Processing of Candidate Object Data 

Additional processing is required before the objects in the 
initial candidate list can be properly analysed. This additional pro- 
cessing is dlagramatlcally Illustrated in Figure 5 and involves both 
geometric correction of the images and navigation (assembly of all 
spacecraft f planet, ephemeris, camera and candidate object information 
into a coherent i^ole) . 

The information obtained in step one (see Figure 2) was used as 
input and was nanually assembled and organised in a file containing the 
image space coordinate of the initial candidate objects and the planet 
center for each relevant approach picture (see Figure 5) , Using a 
file containing the nomirml (l.e., average) locations of the reseau 
markings of the camera, a first order geometric correction was applied 
to all the image space (L,S) coordinates of the initial candidate 
objects and planet centerr. These geometrically corrected coordinates 
arc then called the object space (L,S) coordinates of the initial 
candidate objects and planet centers. 

Using the Viking SEDR files and application programs developed 
for the Voyager spacecraft, the navigation of each image is ecco^- 
lishcd and the orientation Information for the spacecraft, planet, 
camera as well as the object space coordinates of the initial candi- 
date objects is placed in a MIDS (Mars InterMdiate Data Set) cosq>uter 
file. The MIDS file contains all input information necessary for 
further analysis of the pictures. A complete list of all the initial 
candidate objects together with their picture nu^er and their object 
space (L,S) coordinates is given in Table 6 and Table 7 for the Viking 
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1 and Viking 2 approach picturaa, raapactivaly. Additional inforaation 
daacribing the pixel characteriatica, ssixc and brightneaa of the initial 
candidate objecta ia alao included in Table 6 and Table 7 for conven- 
ience. Much of thia additional inforaation vaai however » determined 
later in the analyaia dcacribed below. 

3.4 Analyaia of Candidate Object Data 

The analyaia of the initial candidate objecta list to aearch for 
new aatellitea can be divided into (1) a aiaple global aearch for 
ba.:kground atara and apatial cluatering of candidate objecta and (2) 
a more complex aatellite correlation analyaia of the individual candi- 
date objecta. The global aearch waa performed by properly alipiing 
the planet centcra of all picture framea and properly projecting and 
auperimpoaing all candidate objecta onto one image plane. Thia pro- 
ceaa revealed nothing more than a random distribution of the initial 
candidate object points on the coa^oaite image plane. This implies 
that moat of the initial candidate objecta are random noise and that 
any star in the field of view was too dim to be Imaged. The likely 
magnitude of any star in the field of view should be no more than 
fifth or sixth magnitude which is considerably dimmer than the near 
aeroth magnitude expected for Phoboa (see Table 3) . A more co^lex 
satellite correlation Mialyals* discussed below* waa then developed 
and applied to the Viking 1 and king 2 candidate object data sepa- 
rately. 

The correlation analyaia aearch for satellites is based upon the 
assumption that any real satellite will awve about Mara in accoruance 
with the gravitational inverse-square force law for the planet. The 
object space (L*S) coordinates of those candidate objects occurring 
in different approach pictures that actually correspond to a real 
satellite should therefore coincide when they are dynnically trans- 
formed to a common time frame. A significant effort was expanded to 
develop the application programs suitable to perform this correlation 
analysis. The scope of correlation analysis was restricted to circu- 
r.ar orbits in the equatorial plaiM of Mara. This restriction was 
imposed for two reasons. First, the method could be tested and per- 
fected by locating the satellite Phobo.* idiich has a nearly circular 
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prograde orbit within about one degree of the equatorial plane. Second, 
development and execution of the more general out-of-the-equatorlal- 
plane correlation analysis, although initiated, was soon realized to 
be beyond the time and funding limits of this project. 

Having restricted the satellite correlation analysis to circular 
orbits in the equatorial plane of Mars, the object space (L,S) coor- 
dinates for each candidate object were projected along the optical 
axis of the camera to obtain a unique radius r and longitude 6 for 
that object. Errors in determining the (r, 6) -values exactly arose 
primarily because of the inability to measure the center of Mars 
exactly. The existence of this uncertainty in the (r,d)-value was 
necessarily incorporated in the analysis. The resulting distribution 
in the radii of the circular orbits of the initial candidate objects, 
in 1000 kilometer wide Intervals in the equatorial plane, is summarized 
in Table 8. In Table 8, three of the Viking 2 candidate objects 
(ntimber 40, 41 and 206 in Table 7) have projections that are actually 
inside of the equatorial radius of Mars (3393 km) . 

To minimize the propagation of measurement uncertainties in the 
(r,6) values for each candidate point in the satellite correlation 
analysis, a two-dimensional ideal space (R,0) is Introduced (at a 
reference time t^ midway in the diking approach picture time sequence) 
for searching and correlation display purposes. A computer search 
through this two-dimensional ideal space can then be performed for the 
purpose of correlating the initial candidate object points with points 
in the ideal space. This computer search (l.e., the satellite cor- 
relation analysis) is described in terms of fifteen steps in Table 9. 

The maxlmums of the correlation output arrays X, Y, Z and W defined 
in step 12 provide a direct way of locating the highest correlation 
between initial candidate objects. In step 15, further examination 
of candidate objects producing these correlations allows us to verify 
if any of these correlations correspond to real satellites. 

Vikiftg 1 Correlational Analysis 

The satellite correlation analysis described in Table 9 was first 
applied to the Viking 1 approach pictures, since these pictures are 
significantly mor< ensitlve and less noisy than the Viking 2 approach 
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pictures. Before applying the correlation analysis, 14 initial candi- 
date objects were eliminated from the list. Five of these objects 
[number 2, 7, 42, 49 and 52] were identified as camera blemishes that 
were not in the blemish filter (see Figure 2) and nine objects [number 
17, 18, 19, 20, 21, 32, 36, 38 and 39] were then recognized as strong bit 
errors due to their high DN level above background. The correlation 
analysis was divided into suitable radial regions. For prograde orbits, 
results for the A array are shown in the two-dimensional (rectangular) 
ideal space in Figure 6. The horizontal scale (sample) is the longi- 
tude angle about Mars in 0.5 degree Increments and the vertical scale 
(line) is the radial distance from the center of Mars and ranges from 
8001 to 10,00 j km (from top to bottom) in 3 km Increments. The values 
in the A array are displayed in DN brightness levels, where the bright- 
est possible value is 255 and dark is zero. 

In Figure 6, there are ten long needle-like patterns displayed 
on the ideal space (R,0) plane which actually correspond to the eleven 
initial candidate objects that occur in the radial interval about Mars 
from 8000 to 10,000 km. The transparent overlay identifies each long 
needle-like pattern with the object number of Ita ix.^«;lal candidate 
object given in Table 6 and furthermore shows us that the long needle- 
like patterns for object mutbers 43 and 46 coincide and therefore appear 
to have a single ‘gnature. Note that from Table 8, there are actually 
14 initial candidate objects in Table 6 that are in this radial inter- 
val, but three of these objects (numbers 7, 19 and 32) were previously 
eliminated as blemishes or bit errors (see the discussion above) and 
were thus not include/? fv the correlation analysis. 

Each needle-like pattern is an error ellipse in the ideal space 
for one of the iu^tial candidate objects. Each error ellipse is formed 
by displaying only the 30 top brightness levels of Che A array (l.e., 
values between 255 and 225 DN) . Since this DN level actually repre- 
sents a displacement from the candidate object in the image plane 
(object space) , the error ellipses in Figure 6 represent the possible 
location of each of the 11 initial candidate points in the ideal space 
for a measurement uncertainly of 30 pixels in location of each initial 
candidate object in their approach picture. The major measurement 
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uncertainty is Introduced in determining the planet center and was 
below the 30 pixel error adopted here. Object numbers 3» 5, 10, 14 
and 15, which are properly ordered in event time and all occur before 
the ideal space reference time t^, have negative slopes (as viewed 
from the bottom of Figure 6) . These objects had to be advanced in 
time and the more the objects were advanced in time, the smaller the 
negative slope. Object numbers 27, 34, 43, 46, 48 and 54, which are 
properly ordered in event time and occur after the ideal space refer- 
ence time t^, have positive slopes. These objects had to be backed-up 
in time and the larger the backup time the smaller the slope. Object 
43 and 46 have essentially the same slopes since their event times 
differ only by about 10 seconds (see Table 6 and then Table 1). 

When two or more error ellipses Intersect, this indicates that 
these objects have a high probability of being the same object (l.e., 
a satellite) . Further examination (step 15 of Table 9) is then re- 
quired to more exactly classify the nature of the correlated maximums. 

In Figure 6, object numbers 3, 15, 43 and 46 all Intersect at a common 
point. The only other intersection occurs between object number 46 
and 10. A closeup of these two interaction points is shown in Figure 
7 where brightness contour levels of 245, 235 and 225 DN (i.e., mea- 
surement uncertainties of 10, 20 and 30 pixels) have also been shown. 
Careful examination of the pixel signatures of these initial candidate 
objects shows that object numbers 3, 15, 43 and 46 are several pixels 
in size and that the brightness of the object is consistent with the 
relative camera sensitivity (see Table 6, Table 1 and then Table 5) 
and spacecraft range. Object number 10 is two pixel in size with pixel 
brightness values of 60 and 62 and are hence rather discontinuous with 
the 50 DN background level. This pixel signature is that of a bit error 
and not a real object Imaged by the camera. The intersection of object 
10 and 46 is thus of no significance. 

In Figure 8, the results of the Z-array (see Table 9, step 12) are 
displayed for the intersection point of the object numbers 3, 15, 43 and 
46, and contour levels similar to Figure 7 are also shown. The maximum 
of the Z-array indicates that the object must be located on a circular 
orbit having a radius of 9374 kilometers. This object is the satellite 
Phobos which has an orbital radius between 9376 and 9377 kilometers. 
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A comparison of the four detections of Phobos with other approach pic- 
tures taken In the same time sequences, but with different camera 
parameters. Is given In Table 10. Phobos was detected In addition In 
the picture number 169C65 having a 6-plxel signature with 5 pixels 
having a value of SO DN and 1 pixel having a value of 32 DN above a 
background of 48 DN (a signature not strong enough to trip the object 
locater In Figure 2). In the time sequence 15, the picture numbers 
170C31 and 170C33 had a different field of view (than 170C30 and 
170C32) that excluded Phobos. Transforming In time and space Phobos 
Into all other picture numbers results In either Phobos being out of 
the field of view or too far from the spacecraft to produce a detect- 
able signal. The Viking 1 detection results of Phobos summarized In 
Table 10 also Indicate that all the Viking 2 approach pictures (see 
Table 5) except one [l.e., 219D25 taken In the high gain state that 
effectively doubles the exposure time] will be below the sensitivity 
level required to Image Phobos. 

The correlation analysis was applied to all the Viking 1 Initial 
candidate object data between 4000 and 20,000 kilometer radius. In 
addition to the results for the radial Interval 8000 to 10,000 km dis- 
cussed above, three other correlation maxlmums Involving two Initial 
candidate objects per correlation were also obtained. The first cor- 
relation was between object numbers 9 and 47 having a circular orbit 
of 4967.3 kilometers radius. Object 47 was subsequently shown to be 
a camera blemish and object 7 was not present In the more sensitive 
approach picture 169C47. The second correlation was between object 
numbers 13 and 25 having an orbital radius of 6663 kilometers. Object 
13 and object 23 were not present, respectively. In the more sensitive 
approach pictures (169C47, 169C51) and (170C15) so that the correla- 
tion does not represent a real object. The third correlation was 
between object numbers 8 and 23 having a circular orbit of 10,881 kilo- 
meters radius. Object 23 Is not present In the more sensitive approach 
pictures 170C08 and 170C12 and furthermore. It Is a bit error so that 
the correlation does not represent a real object. 

A retrograde correlational analysis was then applied to all the 
Viking 1 Initial candidate points. Three correlation maxlmums were 
obtained. The first correlation was between object number 9 and 37 


having a circular orbit of 5220 kilometer radius. Both objects were 
subsequently recognized as bit errors. The second correlation was 
between three Phobos objects with object number 3, 43 and 46 giving 
an Incorrect orbital radius of 9156 kilometers. Object 43 and 46 are 
only about 10 seconds apart In event time so that they easily fit In 
a retrograde orbit. The third correlation was between object number 
10 and 27 having an orbital radius of 9528 kilometers. Object 10 Is 
a bit error and object 27 Is not present In the more sensitive time 
sequence approach picture 170C15. The retrograde analysis therefore 
provided no new satellite objects. 

After performing the prograde and retrograde correlational analy- 
sis » the pixel signatures of all the 54 candidates were carefully 
studied and Intercompared . Additional camera blemishes and bit 
errors were Identified and classified. All objects larger than two 
pixels > except Phobos, were Identified as camera blemishes. All re- 
maining two-pixel objects which had a pixel signature larger than 10 
DN above background were seen to be bit errors. In the Viking 1 Ini- 
tial candidate list of 54 objects there are 18 camera blemishes, 19 
bit errors and 4 detractions of Phobos, leaving 13 unidentified two 
pixel objects having brightness signatures no greater than 10 DN above 
background (see Table 6). Because of the Inability to Identify object 
numbers 11, 12 and 14 In the more sensitive time-sequence approach pic- 
ture 160C47, the object numbers 22 and 24 In the more sensitive pic- 
ture 170C08, the object numbers 25 and 27 In the more sensitive pic- 

ture 170C15, object number 33 In the more sensitive picture 170C29, 
object number 50 In the more sensitive picture 170C36 and finally 

object number 51 In the more sensitive picture 170C37, the Viking 1 

Initial candidate list can be reduced to only three unidentified ob- 
jects [numbers 6, 8 and 16]. The two pixel signature of object 6 
(54 and 60 DN above a 50DN background) Is almost certainly a bit error. 
The two-pixel signatures of object numbers 8 and 16 are only 4 DN above 
the background level and are likely random bit errors of low DN values 
but cannot be further differentiated from a satellite object having a 
signal just above the background. 
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viking 2 Correlational Analysis 

Because of the low values of the relative sensitivity and the much 
higher noise levels of the Viking 2 approach pictures In comparison to 
those for the Viking 1 approach pictures and the Inability to detect 
any new satellite In the Viking 1 approach pictures, a complete corre- 
lational an£ lysis for the Viking 2 approach pictures was not performed. 
A limited analysis was, however, applied for prograde orbits to the 
candidate oVects that project Into the equatorial plane In the radial 
Interval between 8,000 and 10,000 km. This radial Interval was chosen 
because It Includes the orbital radius of Phobos (9376-9377 km). Using 
the location of Phobos In the Viking 1 approach picture set, the pre- 
dicted locations of Phobos In all Viking 2 approach pictures were cal- 
culated and compared to location of the Viking 2 Initial candidate 
objects. Only object number 200 In picture number 219D25 compared 
favorably with the predicted location. Picture number 219D23 has the 
highest relative sensitivity of all pictures In the Viking 2 approach 
pictures (see Table 5) because It was taken with a high gain camera 
setting (see Table 2) which doubles the camera sensitivity compared to 
the low gain camera setting. Object number 200 may well be a detection 
of Phobos, but It cannot be confirmed by the remaining less sensitive 
Viking 2 approach pictures. 

The general correlation analysis was then applied to all Viking 2 
Initial candidate objects In the 8000 to 10,000 km radial Interval. 

Five correlations Involving more than two Initial candidate objects 
were present In addition to a number of correlations Involving only 
two Initial candidate objects. Only the five correlations that 
Involved more than two Initial candidate objects were examined. The 
first of these correlations was for an orbital radius of 8040 km and 
Involved object numbers 50, 229 and 247. Objects 229 and 247 do not 
appear In the more sensitive picture number 219D26 of the same time 
sequence so that the correlation does not represent a real object. 

The second correlation, at an orbital radius of 8163 for object numbers 
171, 173, 229 and 247, Is also not real for the same reason since ob- 
ject number 171 does not occur In picture number 219D02 and both object 
numbers 229 and 247 are absent In picture number 219D26. The third 



correlation, at an orbital radius of 8268 for object numbers 30, 38, 

282 and 337, Is also not real. Ojbect number 282 Is a camera blemish 
and the pixel size and brightness of object numbers 337, 38 and 30 
are not consistent when the spacecraft distance from Mars Is considered. 
The fourth correlation, at an orbital radius of 8771 km for object 
numbers 104, 316 and 324, does not represent a real object since ob- 
ject numbers 316 and 324 are not present In the more sensitive picture 
number 219D54. The fifth correlation, at an orbital radius of 8915 km 
for object numbers 141, 183 and 350, Is not of a real object since 
object numbers 183 and 350 are respectively not In the more sensitive 
picture numbers 219D10 and 219D64. 

The final Investigation of Initial candidate objects In Table 7 
was to examine closely all objects larger than two pixels In size. 

This exercise led to the discovery of many new camera blemishes not 
originally contained In the blemish filter (see Figure 2) , to an 
object that was on the planet disc, to an object that was a camera 
reseau, to objects that were bit errors, and finally to objects that 
were random errors (l.e., not visible In a more sensitive picture 
number exposed In the same time sequence). Only object number 278 
which Is three pixels In size, remains as a qualified object larger 
than two pixels In size. This object occurs in picture number 219D44 
which is the most sensitive picture In this time sequence and there- 
fore cannot be disqualified. This object Is likely a random error 
but study of Its pixel signature Is not conclusive because It Is only 
4 DN above the background level. All two-pixel objects In Table 7 
were not analyzed beyond their classification of brightness above 
background . 
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IV. CONCLUDING REKARKS 

Data analysis of the Viking 1 and Viking 2 approach pictures to 
search for new satellites of Mars has been performed using new soft- 
ware tools that were created at the Image Processing Laboratory at 
JPL especially for this project. The search procedure resulted in 
no new satellites. The quality of the Viking 1 approach pictures 
was sufficient to detect objects having similar size and brightness 
as Phobos. Phobos was detected five times in the Viking 1 approach 
pictures and possibly once in the Viking 2 approach pictures. The 
exposure times for the Viking 1 approach pictures were chosen to 
properly image Mars and were significantly shorter than the maxlmxim 
exposure times that would have optimized the camera for detection of 
new satellites. The quality of the Viking 2 approach pictures was 
significantly poorer than the Viking 1 approach pictures because of 
shorter exposure times » inferior transmission of the color filters, 
and noisier background characteristics of the cameras. This made the 
Viking 2 approach pictures essentially useless for detecting new 
satellites of Mars. 

Data analysis of the Viking 1 approach pictures was divided into 
three stages: (1) screening of the image data to obtain an initial 

list of possible candidate objects, (2) navigation processing of the 
candidate object data, and (3) correlation analysis of the candidate 
objects. The first two stages of the data analysis provided a list 
of 54 initial candidate objects (Table 6) which were analyzed in de- 
tail in stage 3. Four of these objects were determined to be Phobos. 

Of the 50 remaining objects, 18 objects were camera blemishes and 19 
objects were bit errors, and 10 objects were identified as random 
errors, leaving only three unidentified objects (object numbers 6, 

8 and 16 of Table 6). Of these three objects, object number 6 is 
almost certainly a bit error and the other two objects are likely 
random bit errors of low DN values which cannot further be differen- 
tiated from a satellite object having a signal just above the background 
brightness level of the image. 
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From the analysis of the Viking 1 and Viking 2 approach pictures. 

It must be concluded that no new satellites of Mars have been discovered. 
This finding Is consistent with the results of Duxbury (1983) ^o, in 
a parallel search effort to discover new satellites of Mars, used pic- 
tures taken by the Viking Orblter late In Its mission. These pictures 
were taken with the Viking Orblter In the equatorial plane of Mars and 
viewed the area within one degree of the equatorial plane and radially 
from the limb of Mars to the orbit of Phobos. Three sequences of pic- 
tures were taken separated by 30 minutes and covered the equatorial 
viewing area on both sides of Mars. Analysis of these pictures by 
Duxbury allowed him to conclude that no satellite with a diameter of 
200 km (or larger) and having a similar albedo to Phobos or Delmos 
were present In his time-space viewing area. The Viking approach 
pictures were much less sensitive but had a superior time window in 
comparison to the Viking Orblter pictures. The two search procedures 
thus had somewhat complementary data sets. 
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Table 1, V iking 1 Approach Picture Information 
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Quality 

View 

A/B 
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Filter 
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DDD 
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407.3 
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VLT 
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135.8 
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08 
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135.8 
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VLT 
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67.88 
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169 
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20 
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1 

169C32 



A 

135.8 
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RED 

563,521 
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169C33 
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407.3 
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VLT 
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23 

17.717 


169C37 



A 

67.88 
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15 

23 
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169C39 



A 

135.8 
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VLT 
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169 

17 

26 

20.760 
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VLT 
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169C46 
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A 
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VLT 

524,555 

169 

17 

26 

47.504 

8 
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A 

407.3 
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RED 
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19 

29 
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67.88 
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169 

19 

29 
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A 
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VLT 
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VLT 
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169C67 
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Table 1 (Continued) 
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170C35 
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170C36 



A 
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VLT 
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170C37 
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57 
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17 
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00 
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B 
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low 

RED 

309,354 

170 

16 

00 

07.536 


170C40 

smear 


A 

407.3 

low 

VLT 

309,342 

170 

16- 

00 

12.152 


170C41 

smear 


B 

135.8 

low 

RED 

309,330 

170 

16 

00 

16.496 

18 

170CA2 



A 

203.6 

low 

VLT 

289,697 

170 

18 

03 

06.235 

i 

170C43 

smear 


B 

135.8 

low 

RED 

289,685 

170 

18 

03 

10.579 


170CA4 

smear 


A 

407.3 

low 

VLT 

289,673 

170 

18 

03 

15.195 


170C45 

smear 


B 

135.8 

low 

RED 

289,661 

170 

18 

03 

19.539 

19 

170C46 



A 

203.6 

low 

VLT 

270,026 

170 

20 

06 

00.318 


170C47 

smear 


B 
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B 
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Table 2 . Vlklns 2 Approach Picture Information 
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Table 2 (Continued) 
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219D37 



A 
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13 

59 

01.384 


219D54 



B 

203.6 

low 

VLT 

223,465 

219 

13 

59 

05.864 


219D57 



A 

33.94 

low 

CRN 

223,430 

219 

7 3 

59 

18.768 


219D58 
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50.91 

low 

CRN 

223,418 
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13 

59 

23.256 


219D59 
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67.88 

low 

RED 

223,405 
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13 

59 

27.745 


219D60 
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lew 

RED 

223,393 

219 

13 

59 

32.225 

18 

219D63 



A 

203.6 

low 

VLT 

203,488 

219 

16 

02 

03.976 


219D64 



B 

203.6 

low 

VLT 

203,476 

219 

16 

02 

08.456 


219D67 



A 

33.94 

low 

CRN 

203,440 

219 

16 

02 

21.811 


219D68 

1 

1 


B 

50.91 

low 

GR.’ 

203,428 

219 

16 

02 

26.299 


219D69 

j 


A 

67.88 

low 

RED 

203,416 

219 

16 

02 

30.788 

i 

j 

219D70 

nissing 


B 

67.88 

low 

RED 

203,404 

219 

16 

02 

35.268 


26 


ORIGIM^L PAdE m 
OF POOR QUALITY 


Table 3. 

Camera Spatial 

Resolution and Phobos 

Imaging 

Spacecraft Distance 
to Mars 
(km) 

Apparent 
Pixel Size 
(km) 

Phobos Image 
Plane Size 
(pixel) 

Phobos 

Brightness 

(magnitudes) 

800,000 

20.0 

1.0 

1.55 

700,000 

17.5 

1.1 

1.26 

600,000 

15.0 

1.3 

0.92 

500,000 

12.5 

1.6 

0.53 

400,000 

10.0 

2.0 

0.42 

300,000 

7.5 

2.7 

-0.58 

200,000 

5.0 

4.0 

-1.46 
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Table 4. Absolute Sensitivity of the Viking 

1 and Viking 

2 Cameras 


Viking 

1 

Viking 2 


Camera A 

Camera B 

Camera A 

Camera B 

Filter Color 

(DN/ms) 

(DN/ms) 

(DN/ms) 

(DN/ms) 

Blue 

9.01 

8.88 

9.98 

7.98 

Minus Blue 

17.84 

14.62 

17.74 

15.35 

Violet 

4.95 

3.82 

3.54 

4.06 

Clear 

26.47 

22.59 

25.07 

22.95 

Green 

13.51 

11.23 

13.24 

11.66 

Red 

5.00 

4.14 

5.06 

4.42 

^These are camera 

filter factors determined for 

the low gain 

state for 


solar irradiance at 1.63AU given in the Viking Flight Memorandum 
0IT-21334-KK (September 21, 1977) 


28 


Approach Pic tures; Summary of the Number of Exposures and their Relative Sensitivity 
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Relative Sensitivity (103 dn) » (Exposure Time) (Filter Factor of Table 4). 


Table 6: Viking 1 i^proach Pictures: Initial Candidate Object List 
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Object Object Blemishes Object Object Brightness above Background (DW) 

Object Picture Coordinates Size and Local Background 

Number Number (L,S) (pixel) Maximums (DN) (1~10) (11-20) (21-30) (.^jQ). 
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Table 7: Viking 2 Approach Pictures; Initial Candidate Object List 
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Object Object Blemishes Object Object Brightness above Background (DN) 

Object Pictures Coordinates Size and Local Background 

Number Number (L.S) (pixel) Maxlmums (DN) (1-10) (11-20) (21-30) (<3Q) 
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Object Object Blemishes Object Object Brightness above Background (DM) 

Object Pictures Coordinates Size and Local Background 

Number Number (L.S) (pixel) Maximums (DN) (1-10) (11-20) (21-30) ( <30) 
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Object Object Blemishes Object Object Brightness above Background (DN) 

Object Pictures Coordinates Size and Local Background 

Number Number (L.S) (pixel) Maxiirums (DN) (1-10) (11-20) (21-30) (<30) 
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Object Object Blemishes Object CBjJect Brightness above Background (DN) 

Object Pictures Coordinates Size and Local Background 

Humber Number (L.S) (pixel) Maximums (DN) (1-10) (11-20) (21-30) (<30) 
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Object Object Blemishes Object Object Brightness above Background (DN) 

Object Pictures Coordinates Size and Local Background 

Number Number (L,S) (pixel) Maximums (DH) (1-10) (11-20) (21-30) (<30) 
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236 219D31 ( 632, 529) 2 - 40-42 

237 219D31 (1042, 397) 2-42 

238 219D31 (1038, 711) 2 - 42-44 

239 219D31 (1049, 914) 2-44 

240 219D31 (1077,1042) 2 - 46-48 


Object Object Bleaishes Object Object Brightness above Background (DN) 

Object Pictures Coordinates Sice and Local Background 

Nunber Nusiber (L.S) (pixel) Max laws (DK) (1-10) (11-20) (21-30) (<30) 
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Object Object Blemishes Object Object Brightness above Background (DN) 

Object Pictures Coordinates Size and Local Background 

Number Number (L,S) (pixel) Maximums (DN) (1-10) (11-20) (21-30) (<30) 
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Table 9. Description of the Satellite Correlation Analysis 


Step Description of Each Step 

1. Choose a reference time to mid-way In the Viking approach 
picture time sequence and define It as the observing time 
appropriate to the Ideal space search analysis. 

2. Choose a particular point (Ro, 0 q) In the ideal space from 
which to Initiate the search procedure. 

3. Project In time and space all Initial candidate object 
coordinates (L,S) Intc their Ideal space coordinates (R, 9). 

4. Select only those candidate objects within ±500 km radial 
distance and ±20^ angular displacement of the Ideal space 
search point (Rq, Oq) . 

5. Time transform the Ideal point (Ro» 9 q) Into the picture 
number frames of each selected candidate object and calcu- 
late the (Lq, Sq) value of the Ideal point in these picture 
number frames. 

6. For each picture number In step 5, compute the distances 
d£ between the Ideal point (Lq, Sq) and each selected 
candidate object (L^, Sj) that occurs In the picture. 

7. Repeat steps 5-6 for all picture numbers containing selected 
candidate objects and build up a set of values [dj^] associ- 
ated with the Ideal space point (Ro» Oq) • 

8. Order the set of values [dj^] and determine the five smallest 
values of d^ [l.e., the five candidate objects nearest the 
point (Ro, 0o)3. 

9. Calculate the quantity (255-dj^) for each of the five values 
(Ignoring negative values) and store them In order of de- 
creasing values in the arrays denoted by A, B, C, D, E, 
where the array elements are those corresponding to the 
location (R^, 0^} In the two-dimensional (rectangular) 
ideal space. 

10. Repeat the steps 2-9 for each point In the Ideal space. 

11. Set a threshold on the arrays B, C, D, and E by deleting 
all Information larger than 30 pixel error In object space 
coordinates (l.e., subtract 225 from the average and Ignore 
all negative values). 

12. Display separately in the two-dimensional (rectangular) 

Ideal space the new arrays [X"A + B, Y-X + C, Z-Y + Z, 
and W * Z + Ej and produce pixel listing and histograms 
(frequency distribution diagrams) of these new arrays. 
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Table 9 (Continued) 


Step Description of Each Step 

13. Use the results of step 12 to manually locate correlated 
maximuras in the two-dimensional (rectangular) ideal space. 

14. Identify the Initial candidate points that produce the corre- 
lated maximums. 

15. Verify these identified candidate points are satellite by 
examining their time and space projection into all picture 
frames and checking that their pixel locations and signatures 
are consistent with the view and camera parameters of each 
picture frame. 

# 


Table 10 


Phoboa In the Viking 1 Approach Pictures 


Phobos Above Spacecraft 


Object 

Number 

Picture 

Number 

Camera 

A/B 

Filter 

Exposure 

(ms) 

Relative 

Sensitivity"^ 

Background 

(DN) 

Range 

(km) 

Satellite 

Imaged 

3 

169C26 

A 

VLT 

407.3 

2.02 

4 

563,591 

Phobos 


C30 

A 

CRN 

67.88 

0.92 

<2 

563,545 

- 


C32 

A 

RED 

135.8 

0.68 

- 

563,521 

- 

15 

169C61 

A 

VLT 

407.3 

2.02 

8 

466,265 

Phobos 


C65 

A 

CRN 

67.88 

0.92 

2-4 

466,218 

Phobos 


C67 

A 

RED 

135.8 

0.68 

<2 

466,195 

- 

43 

170C30 

A 

MBL 

50.91 

0.91 

4 

348,628 

Phobos 

46 

170C32 

A 

VLT 

407.3 

2.02 

8 

348,604 

Phobos 


t 


Relative sensitivity Is defined in Table 5. 
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FIGURE CAPTIONS 
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Figure 1. Viking 1 Approach Picture. A print of the Viking 1 

approach picture number 170C34 is shown that has been 
geometrically corrected and filtered to enhance geolo- 
gical details on Mars. On the planet, circles at 
constant latitudes of 0®, -30", -60® and -85" are indi- 
cated. In the equatorial plane circular orbits with 
radii between 5,000 and 10,000 km are shown at a 1,000 
km spacing. 

Figure 2. Screening of the Image Data. The screening process applied 
to the Viking Orbiter EDR data for the approach pictures 
in order to develop an initial list of satellite candidate 
objects is illustrated. 


Figure 3. Possible Viking 1 Initial Candidate Objects. A print 
of the Viking 1 approach picture 170C32 is shown that 
has been enhanced to show the background brightness levels 
in the image. The locations of possible initial candidate 
objects, numbered both within and beyond the planetary 
disk, are shown in order to facilitate manual screening. 
Object number 347 (upper left) is Phobos. 


Figure 4. Possible Viking 2 Initial Candidate Objects. A print of 
the Viking 2 approach picture 218D24 is shown that has 
been enhanced to show the background brightness levels 
in the image. This image shows greater background varia- 
tion than the Viking 1 picture of Figure 3. The locations 
of possible initial candidate objects, numbered both within 
and beyond the planetary disk, are shown in order to facili- 
tate manual screening. 


Figure 5. First Order Navigation Processing of the Initial Candidate 

Object Data. The procedure applied to geometrically correct 
and navigate the initial candidate object data is illustrated. 


Figure 6. Error Ellipse Diagrams for the Viking 1 Initial Candidate 
Objects. Error ellipses are shown in the ideal coordinate 
frame for the initial candidate objects with radial dis- 
tances of 8001 to 10,000 km (downward vertical axis, 3 km/ 
pixel scale) from Mars and with circular longitude coordi- 
nates from 0 to 360 degrees (horizontal axis, 0.5 degrees/ 
pixel scale) around Mars. The transparent overlay identi- 
fies each long needle-like error ellipse with the object 
number of its initial candidate object given in Table 6. 
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Figure 7. 


Figure 8. 


Error Ellipse Intersections for the Viking 1 Initial Candi- 
date Objects. A closeup of the two Intersections of the 
error ellipses in Figure 6 is shown, where the first three 
overlay brightness contours separated by 10 DN are Indicated. 

Correlation Diagram for Phobos In the Viking 1 Approach 
Pictures. The correlation diagram Is shown for the primary 
intersections of the error ellipses In Figure 6. The 
smallest contour, indicating the most probable location of 
Phobos, provides an orbital radius estimate of 937A kilo- 
meters. See text for discussion. 
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